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Abstract 

A new method for computing hadronic effects on electroweak radiative corrections to low-energy 
weak interaction semileptonic processes is described. It employs high order perturbative QCD re- 
sults originally derived for the Bjorken sum rule along with a large N QCD-motivated interpolating 
function that matches long and short-distance loop contributions. Applying this approach to the 
extraction of the CKM matrix element V u d from superallowed nuclear beta decays reduces the the- 
oretical loop uncertainty by about a factor of 2 and gives V u d = 0.97377(11)(15)(19). Implications 
for CKM unitarity are briefly discussed. 



Precision studies of low-energy semileptonic weak-charged and neutral current processes 
can be used to test the SU(3) C x SU(2)i x U(l)y Standard Model at the quantum loop level 
and probe for potential "new physics" effects. Examples for which a fraction of a percent 
experimental sensitivity has already been achieved include: pion, neutron and nuclear beta 
decays Q], as well as atomic parity violation Q]. In those cases, electroweak radiative 
corrections (RC) have been computed ^, ^[ 0] and found to be significant (of order several 
percent). They must be included in any meaningful confrontation between theory and 
experiment. 

Of course, inherent to any low-energy semileptonic process are uncertainties due to strong 
interactions, since quarks are involved. To minimize such effects, one often focuses on weak 
vector current-induced reactions, where CVC (conserved vector current) protects those am- 
plitudes at tree level from strong interaction corrections in the limit of zero momentum 
transfer. However, even for those amplitudes, electroweak loop corrections can involve weak 
axial- vector effects not protected bv CVC, which give rise to hadronic (strong interaction) 
uncertainties in their evaluation 3j,|4|. In this paper, we focus on the best known and tested 
examples of that phenomenon, the electroweak radiative corrections to neutron and corre- 
spondingly superallowed nuclear beta decays along with their implications for the extraction 
of the CKM matrix element V u d- However, the method we describe is quite general and can 
be easily applied to other charged and neutral current semileptonic low-energy reactions. 

The extraction of V u d (in fact all CKM matrix elements) entails normalizing a semileptonic 
reaction rate with respect to the muon lifetime, or equivalently the Fermi constant derived 
from it 

G„ = 1.16637(1) x 10- 5 GeV~ 2 (1) 

For high precision, electroweak radiative corrections to both processes must be included 
and hadronic as well as environmental effects (e.g., nuclear structure) must be controlled. 
Toward that end, super- allowed (0 + — > + ) nuclear beta decay transitions are very special 
since they only involve the weak vector current at tree level. Small violations of CVC due 
to the up-down mass difference or non-zero momentum transfer are small ~ 0(1CT 5 ) and 
can generally be neglected (or incorporated). Such an analysis leads to the very accurate 
relationship jf],Q| 

lTr ., 2984.48(5)sec /rl , „ , 

1^1 = JIh I D^-n (Superallowed (3 - decays) (2) 
}t{L + KC) 



where ft is the product of a phase space statistical decay rate factor f (which depends on 
the Q value of a specific nuclear beta decay) and its measured half- life t. RC designates the 
total effect of all radiative corrections relative to muon decay as well as QED-induced nuclear 
structure isospin violating effects. It is nucleus dependent, ranging from about +3.1% to 
+3.6% for the nine best-measured super-allowed decays. So, measuring Q and t combined 
with computing RC determines V u d- A similar formula will be given later for neutron beta 
decay. In that case, the Q value = m n — m p is very precisely known, but in addition to the 
neutron lifetime, g& = Ga/Gv must be accurately measured because both weak axial and 
vector currents contribute at tree level jll 0] • 

Our main goal in this paper is to reduce the hadronic uncertainty in the radiative correc- 
tions to super-allowed nuclear beta decays and thereby improve the determination of V u &. 
The need for such an improvement is well illustrated by a survey of ft values and RC for 
super-allowed beta decays by Hardy and Towner more recently updated by G. Savard 
et al. [sj], which found 



where the first uncertainty stems primarily from nuclear structure corrections (including 
0(Z 2 a 3 ) effects) and very small ft value errors while the second, dominant error is due 
to hadronic uncertainties in electroweak loop effects. Although, as we mention later, the 
first error may currently be an underestimate and the central value of V u d could shift due to 
future Q value updates, it is clear that the hadronic loop uncertainty, which comes from weak 
axial-current loop effects, currently limits the determination of V u d and must be improved 
if further progress is to be made. 

Here, we describe a new method for controlling hadronic uncertainties in the radiative 
corrections to neutron and super-allowed nuclear beta decays. It validates our previous 
results jj, Q increasing V u d by only a small +0.00007, but reduces the loop uncertainty by 
about a factor of 2, (0.0004)^^/ — > (0.0002)^^ as we now demonstrate. 

The one-loop electroweak radiative corrections to the neutron (vector current contribu- 
N and snper-aHowednneiearoeta decays are given by Q fl 3 



The first two terms result from loop corrections and bremsstrahlung involving electromag- 
netic and weak vector current interactions, with g(E m ) a universal function [9| integrated 



V ud = 0.9736(2) (4) EW 



(3) 
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over phase space and a short-distance loop effect. They are not affected by strong 

interactions up to ^(^~ - ) — 10~ 5 corrections which can be neglected at our present level 
of accuracy. Higher order leading logs of order a n ln n (mz/m p ) etc. can be summed via a 
renormalization group analysis |4| and 0(Za 2 ) as well as 0(Z 2 a 3 ) contributions have been 
computed for high Z nuclei Q]. They will not be explicitly discussed here, but are included 
in our final results. 

The last three terms in eq. (4) are induced by weak axial-vector current loop effects. 
Their primary source is the 'jW box diagram which involves the time-ordered product of 
the electromagnetic and weak axial- vector currents. That product contains a leading vector 
current component which contributes to + — > + nuclear transition elements. Employing 
the current algebra formulation, one finds |3j 

oo 2 

Box(W)vA= ^J^-^ Fm (5 ) 

where Q is a Euclidean loop momentum integration variable. 

Previous estimates of eq. (5) employed the operator product expansion plus lowest order 
QCD correction to obtain the leading effect Q, Q] 



+ £><£) (6) 



Q 2 [ t 

Integrating over the range m 2 A < Q 2 < oo and combining with smaller vertex corrections 
and ZW box diagrams, that prescription gave a short-distance amplitude contribution 



a 



+ A g 

m A 



, A g ~ -0.34 (7) 



In the numerical estimate, the low energy cutoff was chosen to be m A = 1.2GeV, roughly 
the mass of the Ax resonance, and the error was estimated by allowing to vary up or 
down by a factor of 2. Such a heuristic, albeit crude procedure led to a ±0.0004 uncertainty 
in V U( i- For the long-distance box diagram contribution, nucleon electromagnetic and 



axial-vector dipole form factors were used to find for neutron decay 



as 



C Born (neutron) ~ 0.8g A (/i n + ^ 0.89 (8) 

where Qa — 1-27 and \x n + fi p = 0.88 is the nucleon isoscalar magnetic moment. In the case 
of superallowed nuclear decays, nuclear quenching modifies Cbotu (neutron) and nucleon- 



nucleon electromagnetic effects must be included 



ll|. Overall, in the case of a neutron, 



axial- vector-induced one-loop RC to the decay rate amount to 0.67(8)%. Roughly the same 
uncertainty ±0.08% applies to superallowed nuclear decays. 

To reduce the hadronic uncertainty in RC, we have carried out a new analysis of the 
box diagram axial-vector-induced radiative corrections that incorporates the following 
F(Q 2 ) improvements [12I ]: 

1) Short Distances (1.5GeV) 2 < Q 2 < 00, a domain where QCD corrections remain 
perturbative. 



,2, 1 



W, Q2 



j _ a s{Q 2 )MS _ £<^ a s(Q 2 )MS y _ Q/ a s(Q 2 J]\ls 
71 ~ 71 TX 



(9) 



C 2 = 4.583 - 0.333iV F (10) 
C 3 = 41.440 - 7.607JV F + 0.177JV| (11) 



where Np = number of effective quark flavors. 



2) Intermediate Distances ((0.823GW) 2 < Q 2 < {l.BGeV)' 



Fm = _±12i + e^. . *^ (12) 



m p = 0.776Ge1/ (13) 
m A = 1.230GeV (14) 



1.465GeV (15) 



3) Long Distances: < Q 2 < (0.823Ge1/) 2 

Integrating the long-distance amplitude up to Q 2 = (0.823Ge^) 2 , where the integrand 
matches the interpolating function, and using an update of the nucleon electromagnetic and 
axial-current dipole form factors, we find 

Cbotu (neutron) ~ 0.829 (16) 

a reduction from our own previous result in eq. (8), where the integration was carried up to 
Q 2 = 00. 



Details of the above calculations will be given in a subsequent publication [12j . Here, we 
briefly discuss the results of the above analysis and its implications. 

The QCD corrections to the asymptotic form of F(Q 2 ) have been given in eq. (9) to 
0(a^). The additional terms are identical (in the chiral limit) to QCD corrections to the 
Bjorken sum rule |l3l| for polarized electroproduction and can be read off from well-studied 
calculations [14II15I] for that process. Their validity has been well tested experimentally [16]. 
The interpolating function in eq. (12) is motivated by large N QCD which predicts it 



should correspond to an infinite sum of vector and axial- vector resonances [17( • We impose 
three conditions that determine the residues: i) The integral of eqs. (5) and (12) should 
equal that of eqs. (5) and (9) in the asymptotic domain (l.SGeV^) 2 < Q 2 < oo, which 
amounts to a matching requirement between domains 1 and 2, ii) In the large Q 2 limit, 
the coefficient of the 1/Q 4 term in the expansion of eq. (12) should vanish as required by 
chiral symmetry Q], iii) The interpolator should vanish at Q 2 = as required by chiral 
perturbation theory. Three conditions limit us to three resonances. 

The Q 2 = (0.823GeV) 2 match between domains 2 and 3 was chosen to be the value at 
which eq. (12) equals the integrand of the long-distance contribution. Interestingly, that 
matching occurs near the p mass. A novel technical point in the formulation is that in 
the evaluation of the Feynman diagrams associated with the long-distance contributions the 
integral over the auxiliary variables is carried out first. This leads to integrands that depend 
on Q 2 and can therefore be matched with eq. (12). 

Using this approach, we find that at the one-loop electroweak level the last three terms 
in eq. (4) are effectively replaced by 2.82 - in the case of neutron decay. Comparison with 
eqs. (4) and (8) in conjunction with = 1.2GeV, A g = —0.34 shows that in the new 
formulation these corrections are reduced by 1.4 x 10~ 4 , which increases V u d by 7 x 10~ 5 . 
The smallness of that shift is a validation of our previous result 0, |(| . 

More important than the small reduction in the radiative corrections, our new method 
provides a more systematic estimate of the hadronic uncertainties as well as experimental 
verification of its validity Q|. Allowing for a ±10% uncertainty for the Csom correction in 
eq. (16), a ±100% uncertainty for the interpolator contribution in the (0.823GeV) 2 < Q 2 < 
(1.5GeV) 2 region and ±0.0001 uncertainty from neglected higher order effects, we find the 
total uncertainty in the electroweak radiative corrections is ~ ±0.00038 which leads to ~ a 
±0.00019 uncertainty in V u d- That corresponds to more than a factor of 2 reduction in the 
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loop uncertainty from hadronic effects. 

Employing our new analysis, we find the improved relationship between V u d, the neutron 
lifetime and qa = Gv/Ga 

lTr ., 4908.7(1. 9W , 

\Vud\ = /, , o 2 \ (neutron) (17) 

T n (l + 3g 2 A ) 

Future precision measurements of r n and g& used in conjunction with eq. (17) will ultimately 
be the best way to determine V u d, but for now it is not competitive p. 

In the case of superallowed (0 + — > + transitions) nuclear /3-decays, there are a number 
of corrections, some nucleus dependent, that must be applied to the ft values. They are 
collectively called RC in eq. (2). To make contact with previous studies we factorize 

them as follows: 

1 + RC = (1 + 5 R )(1-5 C )(1 + A) (18) 

The first two factors are nucleus dependent while A is roughly nucleus independent, coming 
primarily from short-distance loop effects. The axial-vector contributions discussed above 



are included in the product (1 + <5r)(1 + A 



orders as well as some next-to-leading logs 



. Because we include leading logs from higher 



IS 



the factorization is not exact and A will 



exhibit some small nucleus dependence. In Table 1, we give the decomposition of RC along 
with their uncertainties, using the results from Hardy and Towner Q, Q for nuclear effects. 
The uncertainty in 1 + 5r comes from Z 2 a 3 and nuclear structure contributions while a 
common ±0.03% error in the Coulomb distortion effect is assigned to 1 — 5c- The entire 
loop uncertainty decribed in this paper is assigned to 1 + A, even though much of it comes 
from medium and long-distance effects and might better be attributed to 1 + 5r. We have 
not made such a distribution, since it is common to all decays. 

Employing the corrections in Table 1 together with eqs. (2) and (18) leads to the V u d 
values illustrated in Table 2. One finds for the weighted average 

V ud = 0.97377(11)(15)(19) (Superallowed (3 - decays) (19) 

Comparing with eq. (3) we see that our analysis gives a somewhat larger V u d due to a 
±0.00007 increase from our new prescription along with refinements from ref. (f| which 
were not included in Savard et al. [8J]. Also, Savard et al. rounded down in their analysis. 
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Table 1. Decomposition of the RC for the nine best-measured supera 



lowed nuclear 
while 1 + 6r 



/3-decays. Coulomb corrections in 1 — 5c are taken directly from ref. 
has been somewhat modified due to our new results. The short- distance 1 + A factor is 
based on the recent update in ref. (f| which includes higher order leading logs and some 
next-to-leading logs. 



Nucleus 


1 + S R 


1 -S c 


1 + A 


10Q 


1.01298(5)(35) 


0.9983(3) 


1.02389(38) 


14Q 


1.01274(8)(50) 


0.9976(3) 


1.02385(38) 


26 Al 


1.01468(21)(20) 


0.9971(3) 


1.02380(38) 


34 C1 


1.01343(34) (15) 


0.9939(3) 


1.02379(38) 


38 K 


1.01322(41)(15) 


0.9939(3) 


1.02378(38) 


42 Sc 


1.01469(49) (20) 


0.9954(3) 


1.02377(38) 


46y 


1.01392(57) (7) 


0.9959(3) 


1.02377(38) 


50 Mn 


1.01394(65)(7) 


0.9957(3) 


1.02376(38) 


54 Co 


1.01398(73) (7) 


0.9947(3) 


1.02376(38) 
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Table 2. Values of V u d implied by various precisely measured superallowed nuclear beta 
decays. The ft values are taken from Savard et al. 8(. Uncertainties in V u d correspond to: 



10 



III, 



1) nuclear structure and Z 2 a 3 uncertainties added in quadrature with the ft error 

n 

2) a common error assigned to nuclear coulomb distortion effects [11], and 3) a common 
uncertainty from quantum loop effects. Only the first error is used to obtain the weighted 
average. 



Nucleus 


ft (sec) 


1+RC 


V ud 






10 C 


3039.5(47) 


1.03542(36) (30) (38) 


0.97381(77) 


;is) 


19) 


14 Q 


3043.3(19) 


1.03441(52)(30)(38) 


0.97368(39) 




19) 


26 Al 


3036.8(11) 


1.03582(30)(30)(38) 


0.97406(23) 




19) 


34 C1 


3050.0(12) 


1.03121(38)(30)(38) 


0.97412(26) 




19) 


38 K 


3051.1(10) 


1.03099(44)(30)(38) 


0.97404(26) 


;i5) 


19) 


42 Sc 


3046.8(12) 


1.03403(54)(30)(38) 


0.97330(32) 




19) 


46y 


3050.7(12) 


1.03376(59)(30)(38) 


0.97280(34) 




19) 


50 Mn 


3045.8(16) 


1.03357(67)(30)(38) 


0.97367(41) 


;i5) 


19) 


54 Co 


3048.4(11) 


1.03257(75)(30)(38) 


0.97373(40) 




19) 






Weighted Average 


0.97377(11) 




19) 
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We note that 46 V gives a somewhat low value for V U d- It differs from the average by 2.7 
sigma. That particular nucleus recently underwent a Q value revision g| which lowered its 
V uc i- It may be indicating problems with other Q values. If the other nuclear Q values follow 
the lead of 46 V, we could see a fairly significant reduction in the weighted average for V u d- 
Clearly, remeasurements of Q values and half-lives of the superallowed decays are highly 
warranted. 

Employing the value of V u d in eq. (19) and the Ki 3 average |3| for V us 

V us = 0.2257(9)(0.961//+(0)), K l3 average (20) 

with / + (0) = 0.961(8) Q leads to the unitarity test 

\V ud \ 2 + \V US \ 2 + \V ub \ 2 = 0.9992(5)^(4)^(8)^(0) (21) 

Good agreement with unitarity is found, with the dominant uncertainty coming now from 
the theory error in the form factor /+(0). Eq. (21) provides an important test of the standard 
model at the quantum loop level and a constraint on new physics beyond the standard model 
at the ±0.09% level. We note, however, that some other calculations j3, Q of /+(0) and 



studies of other strangeness changing decays 



suggest a lower V us value. Combined with 



further Q value revisions possibly leading to a smaller V u d, they could cause a significant 
reduction in eq. (21). A future violation of unitarity is still possible. However, for it to be 
significant, the theoretical uncertainty in /+(0) must be further reduced. 
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